Supersymmetric Phase Transitions and Gravitational Waves at LISA 
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Gravitational waves generated during a first-order electroweak phase transition have a typical fre- 
quency which today falls just within the band of the planned space interferometer LISA. Contrary 
to what happens in the Standard Model, in its supersymmetric extensions the electroweak phase 
transition may be strongly first order, providing a mechanism for generating the observed baryon 
asymmetry in the Universe. We show that during the same transition the production of gravita- 
tional waves can be rather sizable. While the energy density in gravitational waves can reach at 
most Zip !Jgw — 10~ 16 in the Minimal Supersymmetric Standard Model, in the Next-to-Minimal 
Supersymmetric Model, in some parameter range, Iiq f2 gw can be as high as 4 x 10~ n . A stochastic 
background of gravitational waves of this intensity is within the reach of the planned sensitivity of 
LISA. Since in the Standard Model the background of gravitational waves is totally neglegible, its 
detection would also provide a rather unexpected experimental signal of supersymmetry and a tool 
to descriminate among supersymmetric models with different Higgs content. 
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During its evolution the Universe has probably un- 
dergone a series of phase transitions and some impor- 
tant remnants of these events may exist today, includ- 
ing the observed baryon asymmetry of the Universe and 
a stochastic background of gravitational waves (GWs). 
During a first-order phase transition the Universe is 
"trapped" in a metastable state - the false vacuum - 
which is separated from the true vacuum by a barrier 
in the potential of the order parameter, usually a scalar 
field <fi. The transition takes place through the nucle- 
ation of bubbles of the new phase and most of the latent 
heat released in the transition is converted into kinetic 
energy that makes the bubbles expand. When the bub- 
ble walls collide particle production takes place M and 
part of the energy is radiated into gravitational waves 
|pLfl. The properties of this gravitational radiation - 
such as the characteristic frequency and the intensity - 
depend on the typical energy scales involved in the tran- 
sition. In particular, the electroweak phase transition is 
expected to produce a background of GWs with a peak 
frequency around the milliHertz. This frequency happens 
to be the range most relevant for the space interferometer 
LISA 0, which is planned to fly by 2010. 

A strongly first order phase transition could also pro- 
vide a mechanism for generating the observed baryon 
asymmetry in the Universe |tl|. For such a reason, the 
strength of the phase transition has been investigated in 
details in the Standard Model (SM) and in extensions of 
it. Unfortunately, non-perturbative results have revealed 
that there is no hot electroweak phase transition in the 
Standard Model for Higgs masses larger than mw [@- 
Therefore, GWs are not produced at the SM electroweak 
transition. 

Among the possible extensions of the SM at the weak 
scale, its supersymmetric extensions are the best moti- 



vated ones. In the Minimal Supersymmetric Standard 
Model (MSSM), a strong enough phase transition re- 
quires light Higgs and stop eigenstates. For a Higgs mass 
in the range (110 — 115) GeV, there is a window in the 
right-handed stop mass m s t op in the range (105 — 165) 
GeV 13 1. If the Higgs is heavier than about 115 GeV, 
stronger constraints are imposed on the space of super- 
symmetric parameters. However, the strength of the 
transition can be further enhanced in extensions of the 
MSSM, for instance with the addition of a gauge singlet 
in the Higgs sector Q. 

The goal of this paper is to compute the amount 
of gravitational waves generated during the electroweak 
phase transition in supersymmetric extensions of the 
SM. We will find that, depending on the model and 
on the region of parameter space, the stochastic back- 
ground of GWs generated in the collision of bubbles 
nucleated during the supersymmetric electroweak phase 
transition can be within the sensitivity of the planned 
space-interferometer LISA. This opens the exciting pos- 
sibility that the very same supersymmetric physics re- 
sponsible for the generation of the primordial baryon 
asymmetry is also able to provide us with a detectable 
background of gravitational waves. 

A stochastic background of GWs j|,[l5]|ll| can be char- 
acterized by the dimensionless quantity 
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where p gw is the energy density associated to GWs, / is 
their frequency and p c is the present value of the criti- 
cal energy density, p c = 3Hq/(8ttGn), with H Q — 100h 
Km/sec Mpc; ho parametrizes the uncertainty in H . In 
fact, it is more convenient to characterize the stochastic 
background of GWs with the quantity h,Qfl gw (f), which 
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is independent of ho. In Fig. [j]we show the most relevant 
bounds on cosmological stochastic GW backgrounds, to- 
gether with the experimental sensitivities of the vari- 
ous detectors under construction, as discussed in Refs. 
Jl6],[l7]]. In particular, LISA is expected to reach a sen- 
sitivity of order /igfigw — 1CP 12 at / = 1 mHz. At 
this frequency - however - a cosmological signal could 
be masked by an astrophysical background due to unre- 
solved compact white dwarf binaries. Its strength is un- 
certain, since it depends on the rate of white dwarf merg- 
ers and it is estimated to be fo§f2 gw ~ 10~ n [jlO|. At a 
frequency / ~ 10 mHz the LISA sensitivity is expected to 
be of order /igfigw ~ 10~ n , and the astrophysical back- 
ground is below this value. (Correlating two detectors 
one usually gains many orders of magnitudes in the sensi- 
tivity, but because of this extragalactic background, even 
if one would fly two LISA detectors, the sensitivity of the 
correlation would be limited to hgftgw — 5 x 10~ 13 [Q). 
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FIG. 1. The bounds from nucleosynthesis (horizontal 
dashed lines, for N v = 4 and for N v = 3.2), from COBE 
and from ms pulsars, together with the sensitivity of ground 
based detectors and of LISA. See ref. [17] for details. 

Two are the basic quantities which play a role in the 
determination of the GW background generated during 
a first-order phase transition. The parameter a gives a 
measure of the jump in the energy density experienced 
by the order parameter <fi during the transition from the 
false to the true vacuum and it is the ratio between the 
false vacuum energy density and the energy density of the 
radiation at the transition temperature T». The param- 
eter (3 characterizes the bubble nucleation rate per unit 
volume, which can be expressed as T = Toe~xp(/3t) [Q. 
Thus — r/r represents roughly the duration of the 
phase transition and the characteristic frequency of the 
gravitational radiation is expected to be 2wf ~ (3. The 
stronger is the transition, the larger is a and the smaller 
is (3 because of the larger amount of supercooling expe- 
rienced by the system before the transition. The present 
(i.e. properly red-shifted) peak frequency of the GW 
background is determined by [3 and by -ff* , T» , <?* , which 
are respectively the Hubble parameter, the temperature 
and the number of relativistic degrees of freedom at the 
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Typical values for (3/H* for the electroweak phase tran- 
sition are between 10 2 and 10 3 , with T* = 0(100) GeV. 
This gives a frequency / pea k in the range (10~ 4 — 5x 10~ 3 ) 
Hz. From Fig. [j] we see that this is precisely the range 
in which LISA achieves its maximum sensitivity. 

The intensity of the radiation produced is given by H] 
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where v is the velocity at which the bubble expands. In 
the following we will use a value of the velocity v = v(a) 
as given in Ref. Q for bubble detonation. Eq. (||) makes 
it clear that - in order to produce a relevant signal - 
one needs large a and small (3, i.e. a strongly first-order 
transition. 

Before launching ourselves into details, let us briefly 
describe how we have computed the relevant quanti- 
ties appearing in Eqs. (^) and (^|), namely T„, a and 
(3/H*, once the thermal effective potential V((j>, T) for 
the Higgs scalar field(s) is given. The rate of tunneling 
per unit volume from the metastable minimum to the sta- 
ble one is suppressed by the exponential of an effective 

S 3 (T) 

action T = T Q e t , where r is of the order of T , 
and S3 is the extremum of the spatial Euclidean action 

S 3 (T) = J d 3 x [i(V0 h ) 2 +V((p b ,T) computed for the 

configuration of the scalar field(s) describing the bubble 
wall which interpolates between the false and the true 
vacuum. The transition takes place when the proba- 
bility of nucleating a single bubble within one horizon 
volume becomes 0(1). This condition fixes T* to sat- 
isfy S 3 {T*)/T* ~ In (Mpi/100 GeV) 4 ~ 140 for the elec- 
troweak transition. The parameter a is readily computed 
from the definition given above, while [3/H* is given by 



= T, 
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As we have already noticed, extensions of the SM are 
required to obtain a first-order phase transition at the 
electroweak scale. 

In the MSSM two complex Higgs doublets Hi and Hi 
are present in the Higgs sector and new CP-violating 
phases appear which can drive enough amount of baryon 
asymmetry. If the mass tua of the CP-odd field in the 
Higgs sector is much larger than myv , only one light Higgs 
scalar <f> is left and its potential is similar to the one in 
the SM. When tua ~ iri\y the full two-Higgs potential 
should be considered, but the strength of the phase tran- 
sition is weakened (l3| . The one-loop thermal corrections 
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to the effective potential make the quadratic term posi- 
tive at high temperature and create a negative cubic term 
due to loops of the massive bosons in the theory. Due 
to the presence of this cubic term and the positivity of 
the quadratic one, there exists a range of temperature 
in which the point (f> = is a local minimum separated 
from the true simmctry-brcaking one by a small poten- 
tial barrier, that is precisely the set-up for a first-order 
phase transition. The strength of the phase transition is 
enhanced by the presence of new bosons coupled to the 
Higgs, a significant role being played by the right-handed 
stop, which is - apart from the Higgs itself - the lightest 
scalar in the theory and has the largest Yukawa coupling 
to the Higgs 4>. 

To study the amount of gravitational waves gener- 
ated during the electroweak phase transition within the 
MSSM we have made use of the thermal potential cor- 
rected up to two-loop level. Indeed, two-loop corrections 
have been shown to render the phase transition signif- 
icantly stronger in the MSSM |Q. The most relevant 
parameters in the game are the Higgs mass rnhi g g S ) the 
right-handed stop mass TO st0 p and the zero temperature 
ratio between the vacuum expectation values of the two 
neutral Higgses tan/3 M sSM = (H 2 )/(H 1 ). 

Our strategy has been the following. For any given 
choice of the parameters of the model, we have first nu- 
merically computed the nuclcation temperature T„ by im- 
posing that, for the Higgs field configuration describing 
the nucleated bubble, the condition ^(T^/T, ~ 140 is 
satisfied. Then we have computed the parameters a and 
(3 through Eq. (Q). Our results are summarized in Figs. ^ 
and ||. 
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FIG. 2. feo^gw as a function of the stop mass for a 110 
GeV Higgs mass. 

The general prediction is that the intensity of the grav- 
itational waves produced during the MSSM phase tran- 
sition is too small for LISA. For instance, taking a Higgs 
mass of 110 GeV, the right-handed stop mass of 140 
GeV and sin 2 /?mssm = 0.8, we find a~3x 10~ 2 and 
/?/#*- 4xl0 3 , leading to /iQ^gw ^ 2xl0 -18 . Notice that 
one is not allowed to lower too much the right-handed 
stop mass because the thermal squared mass for the stop 
itself would become negative at small temperature, thus 
leading to a physically unacceptable stable color breaking 
vacuum state at zero temperature. 




FIG. 3. fto^gw as a function of the Higgs mass for a 155 
GeV stop mass. 

We have also estimated what happens if we lower the 
Higgs mass down to (the already excluded value of) 80 
GeV, setting the right-handed stop mass at 155 GeV - 
which is the lower value compatible with the absence of 
color breaking minima - and sin 2 (3 = 0.8. We obtain a ~ 
0.1 and /3/H*~ 2xl0 3 , giving /i^ gw ~ 2xl0" 16 , asignal 
still not relevant. The situation does not improve when 
both Higgses are involved in the transition because the 
strength of the phase transition is weaker. A complete 
analysis of the results within the MSSM will be presented 
elsewhere ]l9[ ]. An uncertainty in this estimate is due to 
the determination of v. If the phase transition is not 
strong enough, then v is subsonic, so that the value of 
we have found is really an upper bound. 

The situation improves considerably if we enlarge the 
MSSM sector adding a gauge singlet N Q. This is 
the so-called Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) and is a particularly attractive model 
to explain the observed baryon asymmetry at the elec- 
troweak phase transition. The relevant part of the su- 
perpotential is given by W = XH!H 2 N ~ fiV 3 , where 
now the supersymmetric /i-parameter of the MSSM is 
substituted by the the combination \{N), and k is a 
free parameter. The corresponding Higgs potential reads 
V = Vf + Vd + V so ft, where 

V F = |A| 2 [lAfdffil 2 + \H 2 \ 2 ) + \H.H2l 2 ] 
+ fc 2 |Af| 4 - {Xk*H l H 2 N 2 * + h.c.) , 

V D = ^±Aff 2 | 2 - Iffrl 2 ) 2 + ^\h\h 2 \ 2 , 



V soit = m 2 Hi l^l 2 + m% 2 \H 2 \ 2 + m 2 N \N\ 2 



\AxHxH2N - -kA k N 3 + h.c. 
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The presence of the cubic supersymmetry breaking soft 
terms proportional to the parameters A\ and A k al- 
ready at zero temperature makes it clear that within the 
NMSSM it is quite easy to get a very strong first-order 
phase transition at the electroweak scale [jlif . The or- 
der of the transition is determined by these trilinear soft 
terms rather than by the cubic term appearing in the 
finite temperature one-loop corrections and the preser- 
vation of baryon asymmetry after the phase transition is 
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possible for masses of the lightest scalar up to about 170 
GeV. Barring the possibility that the transition occurs 
along CP-violating directions, the potential becomes a 
function of three real scalar fields (Re TV, Reiii, ReiJ2)- 
In our numerical analysis we have made use of the tree- 
level potential ^) plus the one-loop corrections appearing 
at finite temperature. Overall, we have six free parame- 
ters: the coupling parameters A and fe; the soft-breaking 
mass terms A\ and Ak', the zero-temperature vacuum 
expectation value of the singlet x and tan/^MSSM- 

Given a set of parameters, the strategy has been again 
to determine the nucleation temperature T* and the pa- 
rameters a and (3 which in turn determine the intensity 
and the frequency of the stochastic gravitational back- 
ground. In our analysis we have focussed on those re- 
gions of the parameter space which previous studies have 
shown to give rise to a large enough baryon asymmetry. 



to a rapid growth of the intensity. However, A\, Ak can- 
not be increased beyond the values shown in the figure 
because for higher values the condition 5a(T*)/T* ~ 140 
is never satisfied and the transition does not take place. 

The conclusion is that values /ipfigw — 4 x 10" 11 are 
reachable in the NMSSM in some regions of the parame- 
ter space (we have found similar values of /i^figw also in 
other regions of the parameter space Jll|]). A background 
of this intensity would be within the reach of LISA. 

We find interesting that the very same supersymmet- 
ric phase transition which provides us with a mechanism 
to generate the observed baryon asymmetry might also 
be the source of a sizeable background of gravitational 
waves. If such background is detected, it will be not only 
an indication that supersymmetry might play a role at 
the electroweak phase transition, but even a way to dis- 
criminate among supersymmetric models with different 
Higgs sectors. 
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FIG. 4. hlQ g „ as a function of A\ for A k = 480 GeV, 
x = 350 GeV, A = 0.83, k = 0.67 and tan/3 M ssM = 2 
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FIG. 5. /io^gw as a function of A k for A x = 450 GeV, 
x = 350 GeV, A = 0.83, k = 0.67 and tan/3 M ssM = 2 

Typical results are summarized in Figs. ^ and ^|, where 
we plot /i§i7 gw as a function of A\ and Ak, respectively, 
with A k = 480 GeV in Fig. | and A\ = 450 GeV in 
Fig. [|. Note also that, since the transition in the NMSSM 
is strongly first order, it is correct to use the expression 
for v(a) computed for detonation bubbles |{|. 

The intensity of the produced gravity waves rapidly 
grows as a function of the soft breaking mass terms. This 
is explained by the fact that /ip^gw scales roughly as 
a 3 /(3 2 and a mild increase of a and decrease of leads 
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